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SUMMARY 
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boundary  control  system  of  combined  blowing  and  suction  is  investigated.  Blowing  is  ap- 
plied in  the  front  pa/t  of  the  body  where  the  pressure  gradient  is  favourable,  whereas 
suction  is  applied  in/the  rearigart  of  the  body  where  adverse  pressure  gradients  exist.  In 
order  to  avoid  the  csink  drag"ythe  volume  rate  of  suction  should  be  equal  or  smaller  than 
the  blowing  volume  rate. 

Theoretical  investigations  of  laminar  flows  include  optimization  of  the  blowing  velocity 
distribution  as  well  as  second-order  boundary  layer  effects  such  as  pressure  drag  and  dis- 
placement effects  on  friction  drag. 

Experimental  results  (mean  velocities  and  shear  stresses)  of  turbulent  boundary  layers 
with  very  strong  blowing  velocities  near  the  stagnation  point  of  a circular  cylinder  are 
used  as  a basis  of  a simple  prediction  method.  Experiments  on  a circular  cylinder  show 
considerable  drag  reductions  due  to  the  combined  blowing  and  suction  boundary  layer  con- 
trol system. 


.e  drag 
in  4? 


1.  INTRODUCTION 


The  drag  of  an  immersed  body  can  be  determined  by  integrating  the  forces  acting  on  the 
body  surface.  It  consists  of  two  parts:  the  pressure  drag  due  to  the  pressure  (normal) 

forces  and  the  friction  drag  due  to  shear  (tangential)  forces. 

In  subsonic  flow  the  pressure  drag  is  originated  mainly  by  boundary-layer  separation.  It 
is  well-known  that  boundary-layer  separation  can  be  suppressed  by  boundary-layer  control 
due  to  suction.  The  application  of  continuous  suction  on  the  body  surface  where  usually 
boundary-layer  separation  may  occur  makes  it  possible  to  reduce  the  pressure  drag  of  the 
body. 

On  the  other  hand,  continuous  blowing  perpendicularly  into  the  boundary-layer  leads  to  a 
reduction  of  the  wall  shear  stresses  and  hence  to  a reduction  of  the  friction  drag. 

Figure  1 shows  the  wall  shear  stress  near  the  stagnation  point  of  an  immersed  body  in  in- 
compressible flow  as  a function  of  the  blowing  velocity.  This  result  of  Prandtl's 
boundary-layer  theory  leads  to  the  conclusion,  that  the  wall  shear  stress  can  be  reduced 
to  any  value  wanted  if  the  blowing  rate  is  chosen  strong  enough. 

Reductions  of  both  drag  components,  the  pressure  drag  as  well  as  the  friction  drag,  could 
be  achieved  by  a combination  of  boundary-layer  control  due  to  blowing  and  due  to  suction 
in  suoh  a way  that  blowing  is  applied  on  the  surface  area  with  accelerated  flow  where  high 
wall  shear  stresses  can  be  expected  and  suction  on  surface  areas  with  decelerated  flow 
with  the  danger  of  boundary-layer  separation. 

In  this  paper  theoretical  as  well  as  experimental  investigations  will  be  discussed  which 
refer  to  systems  of  boundary-layer  control  by  combined  blowing  and  suction.  Laminar  as 
well  as  turbulent  flow  will  be  considered. 


2.  LAMINAR  FLOW 

2.1.  Prandtl's  Boundary  Layer  Theory  for  Strong  Blowing  (and  Suction) 

When  the  boundary  layer  alonga  permeable  body  is  controlled  by  blowing  and/or  suction,  the 
solutions  of  Prandtl's  boundary-layer  equations  lead  to  limiting  solutions  for  very  large 
blowing  and/or  suction  rates.  In  Figure  1 the  skin  friction  near  the  stagnation  point  in 
two-dimensional  incompressible  flow  is  shown  as  function  of  the  blowing  velocity  vw 
(negative  vw  corresponds  to  suction) . There  are  two  asymptotes  for  the  skin-friction 
curve  corresponding  to  the  two  limiting  solutions  of  the  boundary  layer  equations.  In  the 
limit  of  strong  suction  the  skin  friction  increases  with  growing  suction  velocity.  The 
skin  friction  is  equivalent  to  the  so-called  "sink  drag"  due  to  momentum  loss  by  suction. 

In  order  to  avoid  this  sink  drag  the  air  sucted  in  has  to  be  blown  out  again  at  a proper 
area.  This  leads  to  another  limiting  solution  for  large  blowing  velocity,  sometimes  called 
the  "blowhard  problem",  see  [1],  As  Figure  1 shows,  in  this  limit  the  skin  friction  is  in- 
versely proportional  to  the  blowing  velocity  vw,  in  other  words,  it  can  be  reduced  as  much 
as  desired,  if  the  blowing  velocity  vw  is  chosen  large  enough.  For  the  given  outer-flow 
velocity  distribution  U(x)  and  the  blowing  velocity  distribution  vw(x)  the  wall  shear  stress 
is  given  by  the  simple  formula: 
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tw(x) 


v U (x) 

Vx) 


dU 

ax 


(1) 


This  formula,  originally  derived  for  blowing  (v  > 0 and  dU/dx  > O)  can  also  be  applied 
for  suction,  if  the  combination  (1/vw(x) ) (dU/dx;  stays  positive,  i.e.  when  suction  is 
used  at  adverse  pressure  gradients  (vw  < 0,  dU/dx  < 0) . As  the  formula  shows  the  wall 
shear  stress  is  inversely  proportional  to  the  blowing  velocity  vw(x).  In  this  limit  of 
strong  blowing  the  flow  field  has  a three-layer  structure  as  shown  in  Figure  2 for  the 
example  of  a circular  cylinder  with  v ^ cos  $,  i.e.  when  blowing  is  applied  on  the  front 
part  of  the  cylinder  and  suction  on  tne  rear  part  of  the  cylinder.  The  viscosity  is  not 
essential  for  satisfying  the  no-slip  condition  at  the  wall  and  plays  a role  only  in  a 
thin  free  shear  layer  along  the  dividing  streamline  detached  from  the  wall. 


The  drag  coefficient  of  the  circular  cylinder  (d:  diameter;  b: 
U„d/v>  is  in  this  limit 

“O'— f 4J- 


span;  R:  radius;  Re  * 
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where 

Vx) 

"57“  /5S  “ CM  008  5 . 


(3) 


The  mass-transfer  coefficient  cM  is  connected  with  the  volume  rate  coefficient 
CQ  " U°b  d 

CD 

by 

cM  “ cQ  /r5 


(4) 

(5) 


where  Q is  the  volume  rate  of  the  recirculating  flow  inside  the  dividing  streamline  (see 
Figure  2) . 

In  Figure  3 the  drag  reduction  of  a circular  cylinder  by  two  different  boundary-layer  con- 
trol systems  is  shown  for  the  Reynolds  number  Re  ■ U d/v  * 4*104.  Curve  I refers  to  the 
combination  of  blowing  and  suction  according  to  Eq.  (3) , whereas  curve  II  corresponds  to 
the  case  of  homogeneous  suction  with  the  total  suction  rate  Q.  The  asymptote  of  curve  I, 
given  by  Eq.  (2) , shows  clearly  that  the  drag  can  become  as  small  as  desired  if  the  cir- 
culating flow  rate  is  large  enough.  On  the  contrary,  curve  II  showB  a drag  reduction  down 
to  a value  of  cD  % 0.2,  but  further  suction  will  increase  the  drag  due  to  the  sink  drag, 
see  also  Figure  1 . The  comparison  shows  very  clearly  that  the  system  of  combined  blowing 
and  suction  is  much  more  effective  than  the  system  of  homogeneous  suction. 

For  bodies  other  than  the  circular  cylinder  the  Eqs.  (2)  and  (5)  will  also  be  valid  with 
just  different  coefficients.  For  a symmetrical  Joukowsky  airfoil  (4.4  % relative  thickness) 
the  drag  coefficient  as  function  of  the  volume  rate  coefficient  is  shown  in  Figure  4.  In 
this  example  the  distribution  of  the  blowing  velocity  v was  chosen  vw  -v  dU/dx  in  order 
to  make  sure  that  the  combination  (1/vw) (dU/dx)  is  always  positive.  As  a consequence,  the 
blowing  rate  Qeiow  is  larger  than  the  suction  rate  Qguc,  namely  QbIow  * 9.1  Qguc.  T^e 
excess  of  blowing  rate  is  equivalent  to  a thrust  of  the  airfoil  similar  to  the  sink  drag 
in  the  case  of  suction.  But  Prandtl's  boundary-layer  theory  is  not  able  to  cover  this 
thrust  effect,  which,  as  will  be  shown  later,  is  a higher-order  boundary -layer  effect. 


2.2.  Optimization 

For  a given  body  geometry  and  a given  velocity  distribution  U(x)  the  following  optimization 
problem  can  be  formulated: 

An  optimal  distribution  vw(x)  has  to  be  found  such  that  the  drag  becomes  a 
minimum  under  the  additional  condition  of  equal  blowing  and  suction  volume 
rates . 

Figure  5 shows  the  result  for  such  an  optimization  of  the  circular  cylinder.  Curves  I 
correspond  to  the  case  shown  in  Figure  2 and  given  by  Eq.  (2).  Due  to  optimizing  the  vw- 
distribution  (curves  II)  the  drag  could  be  reduced  by  44  % from  cD  Cq  Re  = 25  down  to 
Cq  cq  Re  = 14. 

A similar  optimization  of  the  symmetrical  Joukowsky  airfoil  (4.4  % relative  thickness)  led 
to  the  minimum  value  Cp  Cq  Re  = 1.0.  Each  body  has  its  minimum  value  cD  cQ  Re  depending 
on  the  velocity  distribution  U(x). 


2.3.  Second-Order  Boundary -Layer  Theory. 

When  in  the  system  of  combined  blowing  and  suction  the  blowing  rate  is  larger  than  the 
suction  rate,  a negative  drag,  i.e.  thrust,  is  produced.  This  force  against  the  oncoming 
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flow  direction  must  be  a pressure  force.  This,  however,  cannot  be  predicted  by  Prandtl's 
boundary  layer  theory.  This  thrust  due  to  blowing  is  a second-order  boundary  layer  effect. 
Second-order  boundary  layer  theory  will  also  yield  information  about  the  displacement 
effects  of  the  fairly  thick  boundary  layers  due  to  blowing. 


The  second-order  boundary-layer  effects  for  the  circular  cylinder  (vw  <\>  cos  <j>,  see  Figure 
2)  have  been  investigated  in  [2].  The  drag  formula,  Eq.  (2),  has  been  extended  to 


= 8tt  + 2n2 
cM/Re  Re 


(6) 


It  is  worth  mentioning,  that  the  second-order  term  is  independent  of  cM. 

Figure  6 shows  the  second-order  effects  on  the  drag  of  the  Joukowsky  airfoil  (vw  ^ dU/dx, 
see  Figure  4)  at  Reynolds  number  Re  =>  U^c/v  “ 104  (c  ■ chord  length  of  airfoil).  The 
blowing  rate  cqbiow  is  larger  than  the  suction  rate  cQSuJ,  therefore  thrust  due  to  blowing 
occurs,  cD  = - 1.78  cqbiow  As  example  of  the  circular  cylinder  an  additional  fric- 

tion drag  due  to  second-order  effects  {displacement  effects)  exists, which  reaches  a finite 
value  of  cD  % 0.0052  (independent  of  cqbiow)  for  infinitely  large  blowing  rates.  Contrary 
to  the  circular-cylinder  example  where  second-order  effects  on  drag  were  only  friction 
drag,  here  also  a pressure  drag  (curve  lib)  exists,  which  unfortunately  has  the  same  order 
of  magnitude  as  the  friction-drag  reduction  due  to  blowing. 


3.  TURBULENT  FLOW 


3.1.  Experiments  of  the  Flow  Field  Near  Stagnation  Point 


So  far,  laminar  flow  was  assumed.  But  in  the  area  where  blowing  is  applied,  turbulent  flow 
has  to  be  expected,  because  blowing  has  a strong  destabilizing  effect  on  the  boundary  layer. 
Therefore,  detailed  measurements  of  the  flow  characteristics  (mean  velocities,  turbulence 
intensities,  turbulent  shear  stresses)  within  the  boundary  layer  with  blowing  have  been 
undertaken  by  using  Laser-Doppler  anemometry,  [4].  The  experiments  have  been  carried  out 
at  the  stagnation  point  area  of  a porous  circular  cylinder  (see  Figure  2)  with  extremely 
large  blowing  rates  (up  to  vw/U00  * 0.18)  which  have  not  yet  been  reported  in  the  literature. 

Figures  7 and  8 show  the  mean  velocities  and  the  turbulent  shear  stresses  at  different  po- 
sitions over  the  front  part  of  the  circular  cylinder.  Turbulent  boundary  layers  along  im- 
permeable walls  have  usually  a two-layer  structure  with  the  "wall"  layer  and  the  "defect" 
layer,  where  the  friction  velocity  ux  = Ay/p  is  the  characteristic  velocity. 


Turbulent  boundary  layers  with  strong  blowing  still  have  the  two-layer  structure,  but  with 
a few  important  differences. 

Firstly,  the  basic  characteristic  velocity  is 


u _ V Tmax  ' 

t max  p 


(7) 


i.e.  based  on  the  maximum  shear  stress  t max  within  the  boundary  layer  (see  Figure  8) 
rather  than  based  on  the  wall  stress  tw. 

Secondly,  the  boundary  layer  thickness  is  growing  with  the  blowing  velocity,  but  the  de- 
fect layer  is  growing  relatively  faster  than  the  wall  layer,  in  other  words,  the  defect 
layer  plays  an  increasing  role  compared  to  the  wall  layer  when  the  blowing  rate  is  in- 
creasing. 

As  a consequence,  all  velocity  distributions  shown  in  Figure  7,  will  approximately  collapse 
into  one  universal  curve  (except  close  to  the  wall)  if  represented  as  defect  laws  according 
to 


U - u 
ut  max 


-•  <i> 


(8) 


see  for  details  [4] . This  law  can  be  interpreted  as  an  extension  of  the  law  of  the  wake  by 
McQuaid  [5].  The  function  g (y/6)  is  identical  with  the  one  found  by  McQuaid  for  blowing 
rates  up  to  vw/U  = 0.008.  It  should  be  mentioned  that  McQuaid  used  a characterstic  velocity 
different  from  uT  max.  In  [4]  however  it  is  shown  by  asymptotic  arguments  why  uT  max  is 
the  proper  velocity  scale. 


In  Figure  9 the  experimental  results  are  shown  as  defect  laws  according  to  Eq.  (8)  and  com- 
pared with  McQuaid' s defect  law.  Particularly  the  results  for  <f>  ■ 30°  show  that  the  de- 
fect law  covers  most  of  the  boundary  layer  for  cases  of  strong  blowing. 


3.2.  Prediction  of  Turbulent  Boundary  Layers  with  Strong  Blowing. 

Turbulent  boundary  layers  with  strong  blowing  are  very  well  represented  by  th-  defect  layer 
governed  by  Eq.  (8)  . This  result  was  used  as  a basis  for  a simple  prediction  method  of 
such  boundary  layers.  The  following  equations  have  been  applied: 


1-M 


a)  Momentum  integral  equation  for  the  defect  layer  (neglecting  the  small  wall  shear  stress): 

~(9) 


S1  » S + S l52  °2>  - »»  0 • 


b)  Velocity  distribution  according  to  defect  law  of  Eq.  (8) : 
u 


u _ . t max 
U U 


g ( 6“) 


(10) 


This  leads  to 


51  _ f ,ut  maX' 
6 U ' 


u 

2 _ * , t max, 

T “ f2l“u 


(ID 


There  are  two  unknowns  6(x)  and  uT  max^x^  for  given  U(x)  and  vw(x).  In  addition  to  Eq.  (9) 
an  entrainment  relationship  has  been  used: 


dQ  _ d_ 
dx  dx 


[U  (6  - 6^] 


VE  + VW 


(12) 


where 


t max 


= c„  <■  0.08 
E 


(13) 


a value  typical  of  entrainments  of  wakes. 

The  Eqs.  (9)  to  (13)  are  a complete  set  to  predict  6 (x)  and  uTmax(x)  forgiven  U(x)  and 
vw(x).  Figure  10  shows  the  results  of  such  a prediction  and  the  comparison  with  the  ex** 
periments  described  in  Chapter  3.1.  The  wall  shear  stress  depicted  in  Figure  10c  has 
been  determined  from  uT  max(x)  by  a simple  estimation  given  in  [4]. 

The  prediction  method  presented  here  can  be  used  for  a good  estimation  of  the  behaviour 
of  turbulent  boundary  layers  with  very  strong  blowing. 


4.  EXPERIMENTS  ON  DRAG  REDUCTION  OF  A CIRCULAR  CYLINDER 

The  combination  of  boundary  layer  control  by  blowing  and  suction  has  been  applied  to  a 
circular  cylinder  in  incompressible  flow.  In  Figure  11a  typical  arrangement  is  shown.  The 
front  quarter  of  the  cylinder  periphery  was  used  for  blowing,  whereas  the  rest  of  the  peri- 
phery was  used  for  suction.  The  suction  velocity  was  one  third  of  the  blowing  velocity  to 
make  sure  that  the  net  volume  rate  was  equal  to  zero.  Typical  results  for  the  pressure 
distributions  at  Re  = 9 * 104  are  shown  in  Figure  12.  Due  to  boundary-layer  control  the 
pressure  drag  could  be  reduced  considerably  from  cp  = 1.04  without  control  down  to  cD  = 0.52 
for  vw  sue  /U®  = 0.03  . From  the  comparison  with  the  pressure  distribution  of  potential 
theory  it  can  be  concluded  that  even  better  results  might  be  possible  by  further  refine- 
ment of  the  combined  boundary  layer  control  system. 
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Figure  1 : Skin  Friction  near  the  Stagnation  Point 

in  Two-dimensional  Incompressible  Flow 


Figure  2:  Flow  past  a Circular  Cylinder 

Blowing  in  the  Front  Part 
Suction  in  the  Rearward  part 

Three  Layer  Flow 

I : Inviscid  Outer  Flow  (potential  Flow) 

II  : Free  Shear  Layer 

III:  Inviscid  Inner  Layer  with  Vorticity 
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Figure  3:  Drag  Reduction  of  a Circular  Cylinder  by  Two 

Different  Boundary-Layer  Control  Systems. 
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I : Combined  Blowing  and  Suction  vw  - cos  $ , 

cn  : Volume  Rate  Coefficient  of  Recirculating 

U _4 

Flow.  Asymptote:  cD  = 6.25  • 10  /c^ 


II:  Homogeneous  Suction,  vw  = const. 
Cq  : Suction  Coefficient 
Asymptote:  c„  = 2 • cn  (Sink  Drag) 
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Finite  Difference  Solution 


v ^-v  ^Asymptote 
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^cD  = -2(c08|ow-lc0JI=-1.7Bca#)e. 


Figure  6:  Drag  Reduction  of  Joukowsky  Airfoil 

by  Combined  Blowing  and  Suction. 


I : Prandtl's  (first  order)  Boundary-Layer  Theory. 

-4 

Asymptote  for  Friction  Drag:  cQ  * 9.1*10  /cQ  b1qw 
cq  Blow  : B1°win9  Rate  Coefficient  =9.1  cQ  guc 


II:  Second  Order  Boundary-Layer  Theory. 
Ila:  Friction  Drag 
lib:  Pressure  Drag 

Thrust  due  to  Blowing,  cp 
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Figure  7:  Mean  Velocities  in  the  Boundary  Layer  in  the 

Stagnation  Region  of  a Circular  Cylinder  with 
Strong  Blowing  (LDA  Measurements)  Re  = Uwd/v  , 
U * U(x)  : Local  Outer  Flow  Velocity. 
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